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Abstract 

We discuss the relationship between exact solvability of the Schròdinger 
equation with a position-dependent mass and the ordering ambiguity in 
the Hamiltonian operator within the frame of supersymmetric quantum 
Q_i| mechanics. The one-dimensional Schròdinger equation, derived from the 

general form of the effective mass Hamiltonian, is solved exactly for a 
system with exponentially changing mass in the presence of a potential 



fi 



with similar behaviour, and the corresponding supersymmetric partner 



Hamiltonians are related to the effective-mass Hamiltonians proposed in 
the literature. 



"íh ■ PACS Numbers: 03.65Ca, 03.65Fd, 03.65Ge 



1 Introduction 

The study of quantum systems with position-dependent effective masses has been 
the subject of much activity in recent years. The Schròdinger equation with non- 
constant mass provides an interesting and useful model for the description of 
many physical problems. The most extensive use of such an equation is in the 
physics of semiconductor nanostructures. This field has arisen due to the im- 
pressive devolepment of sophisticated technologies of semiconductor growth, like 
molecular beam epitaxy, which made it possible to grow ultrathin semiconductor 
structures, with very prominent quantum effects 0. The motion of electrons in 
them may often be described by the envelope function effective-mass Schròdinger 
equation, where the material composition- (i.e., the position-) dependent effec- 
tive mass of electrons replaces the constant particle mass in the conventional 
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Schròdinger equation. The most popular of these structures is the semiconductor 
quantum well, and the Schròdinger equation here is effectively one-dimensional. 
Another instance where such an equation is employed, this time three-dimensional 
and with spherical symmetry, is in the pseudopotential-theory-based density func- 
tional calculations in sòlids: to reduce the còmput ational load, model pseudopo- 
tentials with position-dependent electron mass which replace nonlocal pseudopo- 
tentials have been considered 0. 

Since the momentum and the mass operators no longer commute in case of 
spatially varying mass, a question concerning the correct form of the kinetic en- 
ergy operator of the generalized Hamiltonian has arisen. This problem of ordering 
ambiguity is a long standing one in quantum mechanics, see for instance the ex- 
cel·lent critical review by Shewell pi. There are many examples of physically 
important systems, for which such ambiguity is quite relevant. For instance we 
can cite the problem of impurities in crystals [|, [5|, |B| , the dependence of nuclear 
forces on the relative velocity of the two nucleons |7|, ||, the minimal coupling 
problem in systems of charged partides interacting with magnètic fields ||, and 
more recently the study of semiconductor heterostructures |], [10], 11 1. 



Notwithstanding, taking into account the spatial variation of the semiconduc- 
tor type, some effective Hamiltonians are proposed with a spatially dependent 
mass for the carrier [12-17], and many authors have been trying to determine 
the correct Hamiltonian phenomenologically. In this article we try to circum- 
vent the problem of ambiguity by presenting a scheme to obtain unambiguously 
the Schròdinger equation with spatially varying particle mass, which makes clear 
the link between possible choices of the kinetic energy operator for quantum 
systems with position-dependent effective mass, within the frame of supersym- 
metric quantum mechanics |[18|| . The strategy here is to tackle the problem with 
a fundamental point of view, ie., without using a particular form of the effec- 
tive potential. Within this framework and using the supersymmetric formalism 
we will show that one can arrive at a conceptually consistent result for exactly 
solvable systems. 

The application of supersymmetry ideas to nonrelativistic quantum mechanics 
has provided a deeper understanding of analytically solvable Hamiltonians, as well 
as a set of powerful approximate schemes for dealing with problems admitting no 
exact solutions. The concept of shape invariance |Tí| has played a fundamental 
role in these developments. The aims of the present work are to consider the 
application of the supersymmetric approach to quantum systems with position- 
dependent mass and to extend the concept of shape invariant potentials to the 
nonconstant mass scenario to see clearly the relation between the effective-mass 
potentials existing in the literature and supersymmetric partner potentials. 

The paper is organized as follows. In section 2 we provide a brief review of the 
Schròdinger equation for systems with position-dependent effective mass. Section 
3 deals with the exact solvabilty of a system involving a particle with a spatially 
dependent mass in an arbitrary potential well. The application of supersymmetric 
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approach to this system, together with the use of shape invariance concept to 
deduce the correct operator ordering for the Hamiltonian, are also studied in this 
section. Finally, some conclusions are drawn in section 4. 



2 Generalized Schròdinger equation 

We start this section by defining a quite general Hermitian effective Hamiltonian 
for the case of a spatially varying mass which will be denoted by m [= m(x)} 
throughout the present work. In general, the Hamiltonian proposed by von Roos 



20|j is used, 

HvR= 4 



m a pmP p m? + m? p vnP p m a 



(1) 



but to accomodate the possibility of including the case of the Weyl ordering [fIT] 
in a more evident way, we will use an effective Hamiltonian with four terms given 
by @ 

H = — — ^ |a [m~ 1 p 2 + p 2 m~ 1 ] + m a pmP pm? + m? pmP 'pm a } . (2) 

In both cases, Eqs. (1) and (2), there is a constraint over the parameters such 
that a + (3 + 7 = —1. 

The one-dimensional time-independent Schròdinger equation reads 

(h 2 \d 2 ^ Ti 2 (m'\dïjj ÍTT 1 

where U aia involving all the ambiguity is 

fi 2 r , 

(a + 7 — a)m m" + 2(a — «7 — a — 7)777. , (4) 



u <*yaW 4m 3 (a + 1) . 

in which the first and second derivatives of m(a;)with respect to x are denoted by 
m'and m". It is clear that the effective potential is the sum of the real potential 
profile V(x) and the modification U aia emerged from the location dependence 
of the effective mass. A different Hamiltonian leads to a different modification 
term, see Table 1. 

It is curious to note that all the ambiguity is in the U aia term, and that it 
can be eliminated by imposing some convenient constraints over the ambiguity 
parameters, namely 

a + 7 — a = 0, a — «7 — a — 7 = 0, (5) 

which have two equivalent solutions,(i) a = and a = 7, or (ii) a = a and 
7 = 0. In this case the effective Schròdinger equation will not depend on the 
ambiguity parameters, but will contain a first order derivative term. In the next 
section, we will be interested in getting exact solutions of the resulting equation 
for a particular potential, and trying to get some information about the proposed 
orderings appearing in the literature. 
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3 An exactly solvable system 



The interest in exactly solvable problems in quantum physics has increased sharply 
in the last few years. This is concerned, of course, with the fact that the de- 
scription of the behaviour of nonconservative physical systems is usually very 
complicated, but in some cases such systems can be modelled by means of quite 
a simple Hamiltonian, which leads to Standard problems of quantum mechanics. 

Starting with the Schròdinger equation in Eq. (3) and making a new defmition 
for the wave function 

ip(x) = m l l 2 (ç(x), (6) 
one gets a differential equation in a more familiar form 



h 2 d 2 <p 



+ {U eff -E)<p=0, 



2m dx 2 

with a new effective potential defined through 

h 2 

U eff{x) = V(x) + U a7a (x) + — 



(7) 




m 
m 



(8) 



To demonstrate the simplicity of the present method, we consider here a partic- 
ular case, which has been recently studied p2| , where one have exact solution for 
the above equation. That is a particle with exponentially decaying or increasing 
mass in the presence of a potential with similar behaviour, 



m(x) = m e cx , V(x) = V e c 



(9) 



where rriç, is a constant mass. This problem is often encountered in the calculation 
of confmed energy states for carriers in semiconductor quantum well structures 
under the envelope-function and the effective-mass approximations where the ef- 
fective mass of a carrier is spatially dependent on the graded composition of the 
semiconductor alloys used in the barrier and the well region of the microstruc- 
tures. 

Multiplying each term in Eq. (7) by m(x) and dividing by tuq, one arrives at 
an usual Schròdinger equation for the system of interest, 



h 2 d 2 <p 
2m dx 2 



+ (Vc 



o e 



2cx 



where 



h 2 

— (1-c 2 / 
m 



8), q 



E e cx ) <f = e<f 



■(a — 2a / y — a — 7). 



(10) 



(11) 



4(a + l) 

Note that Eq. (10) corresponds to a Schròdinger equation for a particle with con- 



stant mass under the influence of the Morse potential PB[. It is thus clear that 
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if one knows the spectral properties of the constant-mass Schròdinger equation 
of any potential, one then can readily obtain a corresponding potential for the 
effective-mass Schròdinger equation with identical spectral properties. Proceed- 
ing with the well known energy spectrum of the Morse potential, 



2„2 



2mn 



VïmoVo ( , 1 
— [n + - 

hc V 2 




Q- 



(12) 



from which we find the energy spectrum of the effective potential appearing in 
Eq. (7), that is 



E,, 



h c 



2mr 



[2n + 1 + i/(a,7, a)]. 



where the ordering term v{a, 7, a) is 



i/(a,7,a) 



a + 1 



2«7 



a 



7). 



(13) 



(14) 



One can now study the effect of using some of the orderings appearing in 
the literature. Considering Table 1, it is not difficult to see that the ambiguous 
term v is zero for the effective Hamiltonians in Refs. [15, 16, 
have different orderings, while v — 1 for the BenDaniel-Duke Hamiltonian 
For both cases {y = 0, 1), the corresponding Hamiltonians have exactly the same 

nH have one 



although they 



15, lí 



spectra except for the fact that the Hamiltonians in Refs. 
bound state more than the BenDaniel-Duke Hamiltonian. Thus, they can be 
treated as the supersymmetric partner Hamiltonians which is the subject of 
the next section. Furthermore, one ends with a complex energy for the orderings 
proposed in the literature |T~ 
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which could be possibly discarded due to 
the physically unacceptable energies. This makes clear that unacceptable physics 
consequences may occur unless specific choices are made in the Hamilton operator 
ordering for a system undertaken. 

As we are dealing with a confmed particle system, one may wish also to con- 
firm Eq. (13) by mapping the Morse potential onto harmònic oscillator system, 
which seems more realible than the Morse oscillator for the system of interest. 
For this purpose, we invoke the change in the variable as well as in the wave 
function, 

W lc ^(x) = J2fcy- l / 2 F{y) . (15) 



x 



This reduces Eq. (10) to an equivalent Schròdinger equation 



k 2 d 2 F 
2m dy 2 



+ 



W 2 



h z c 2 + 32m e\ 
8m c 2 !/ 2 ) 



AE 

— F ' 



(16) 



in a more familiar form involving a harmònic oscillator potential with centripetal 
barrier. In contrast to Eq. (10), which contains a variable parameter e repre- 
senting the Morse oscillator energy, we have a constant term on the right-hand 
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side of Eq. (16), and the energy term (e) is contained in the effective potential 
parameter. Thus Eq. (16) may be looked upon as the radial Schròdinger equation 
with a fixed energy but with a variational angular momentum quantum number. 
From which one can easily arrive at Eq. (13), which clarifies that both treatments 
(Morse and harmònic oscillator mapping) are equivalent to each other. 

In the following section, we will focus our attention on how to apply the 
supersymmetric quantum mechanical formalism to the system under considera- 
tion in order to clear out the hidden relation between the effective Hamiltonians 
proposed in the literature for the spatially dependent mass. 



3.1 supersymmetric approach 



The problem of generating isospectral potentials in quantum mechanics has been 
considered for more than 50 years, but recently the research efforts on this tòpic 
have been considerably intensified. A new field, supersymmetric quantum me- 
chanics (SUSYQM), devoted to this class of problems has emerged, which deals 
with pairs of Hamiltonians that have the same energy spectra, but different eigen- 
states. A number of such pairs of Hamiltonians share an integrability condition 
called shape invariance JÏ9[ . Although not all exactly solvable problems are shape 



invariant [^4]], shape invariance, especially in its algebraic formulation, is a pow- 
erful technique to study exactly integrable systems which have always been at 
the centre of attention in physics and mathematics. 

It would be interesting therefore to extend the SUSYQM to handle cases with 
position-dependent mass. Recently supersymmetric techniques have been applied 
to obtain exact solutions of Schròdinger equation with nonconstant mass |25], ^(J . 
Using the spirit of these works, here we generalize the supersymmetric formalism 
for the problem considered in the previous section. All considerations are made 
for the one-dimensional Schròdinger equation. 

Proceeding as in the case of constant mass, we introduce a superpotential 
W(x) and the associated pair of operators A and A + defined by 



Aip 



h dip 
/ 2m dx 



+ Wip 



A + = ~ 




+ Wip 



(17) 



Notice that, due to the position dependence of the mass, d/dx and %/ \j2m do 
not commute anymore. Within the framework of SUSYQM, the first partner 
Hamiltonian reads 

„, = A+A = _«L* . f f]' * + \ w i _ f «L v 

2m dx 2 \ 2m J dx y \ v 2m ) 

where the prime denotes the first derivative with respect to the variable x. 
At this stage, we propose an ansatz for the superpotential, 

% C r- fi C 



W(x) 



2m 



m 



(19) 
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H = ^+ = -— — -f— Y— + 

2m dx 2 \ 2m J dx 



and note that the Hamiltonian in Eq. (18) found via the supersymmetric formal- 
ism corresponds to the specific effective-mass Hamiltonians in Eq.(3) with real 
potential profile V(x). Explicitly, these are the Zhu and Kroemer Hamiltonian 
(a = 0, a — 7 = — 1/2) [jnj, the Li and Kuhn Hamiltonian (a = a = 0,/3 = 7 = 
— 1/2) ||16|| , and the Weyl Hamiltonian (a = 1, a — 7 = 0) pi| , in which the effec- 
tive potentials (for which u(a, 7, a) = 0) may be expressed in the supersymmetric 
form, 

V eff (x)=V e™ + U aja (x) = W 2 -^j . (20) 

This justifies the operator ordering in the Hamiltonian used by many authors 
in the rough calculation of the confinement states in quantum well structures 
in the effective mass scheme, as the first partner leads to the exactly solvable 
Hamiltonian systems within the frame of SUSYQM. 

The associated supersymmetric partner Hamiltonian, for the case v(a, 7, a) = 
1 which describes the BenDaniel-Duke effective mass Hamiltonian (a = a = 7 = 
0) 0, is 

w2 _ í nw y | 2nw _ í h y 

(21)" 

We see that the two supersymmetric partner Hamiltonians Hi and H 2 describe 
partides with the same effective mass-spatial dependence, but in different poten- 
tials. The second partner potential corresponding the BenDaniel-Duke effective 
potential does not incorporate the ambiguity term, 

Vbdd(x) = y <■ = W< _ (W)' f » V , ,22) 



'2m j 

where the double prime denotes the second derivative with respect to x. 

There is a correspondence between the energy eigenvalues of the isospectral 
Hamiltonians H\ and H 2 , although they have different effective potentials. The 
energy of the nth bound state of Hi coincides with that of the (n — l)th bound 
state of H 2 , which is the case expressed through Eq. (13). The ground state 
of H% has no associated state of H 2 . The reader is referred to Ref. [Ï8[] for a 
comprehensive review of the supersymmetric quantum mechanics. 

For an exactly solvable system it is thus obvious that a proper choice of the 
effective mass and the corresponding superpotential would impose the elimination 
of the supersymmetric partner potential that involves all the ambiguity. This 
makes clear that which choice of the effective mass Hamltonians proposed in the 
literature is physically acceptable for the system undertaken. 

Finally we note that, from the relation between the superpotential and the 
ground state wave function of Hi, 

W {X) = (23) 
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it is straightforward to solve ipi =0 , which satisfies the required boundary condi- 
tions dictated from the conservation of current through the envelope functions of 
heterostructures. 

In sum, considering the definition of the effective potential which is the sum- 
mation of a real potential V(x) and a term U aia (x) resulting from the mass 
dependence on the location, and bearing in mind that the effective potential re- 
fies on the Hamiltonian utilized, we have connected the deviation from the real 
potential, due to the ordering ambiguity, to the supersymmetric partner Hamil- 
tonians. It is found that the operator ordering in the kinetic energy operator to 
be the same as that endorsed by recent authors [27, and the references therein] 
for the calculation of confined states in semiconductor microstructures under the 
simplified effective-mass and envelope- function, with a correctional term U aia (x) 
that is proportional to the derivatives of the mass profile m(x). 



4 Concluding remarks 

In this work we have first discussed the problem of solvability and ordering am- 
biguity in quantum mechanics, as the form of the effective mass Hamiltonian has 
been a controversial subject due to the location dependence of the effective mass. 

It was shown through a particular example that exact solutions could be used 
as a kind of guide, at least, restrict the possible choices of ordering. The principal 
idea is to suppose that once one have found the ordering without ambiguity for 
a given potential or class of potentials, that ordering should be extended to 
remaining physical potentials. 

Our work has also made clear that the Hamiltonian proposed by Li and Kuhn 
[ y| is in fact equivalent to that coming from the Weyl ordering pï| , as can 



be easily checked from Eqs.(4, 21, 24) the effective potentials and consequently 
from Eqs. (13,14) the energy spectra of both these orderings are equal. In 
addition, though the ordering in the Zho and Kroemer effective Hamiltonian 



operator |Ï5| is different than those of Li and Khun, and that of Weyl, we have 
arrived at the same result in each case within the frame of both, the Standard 
quantum theory and supersymmetric quantum mechanics. Hence, we remark 
that when the ambiguity term is linear in the momentum, contrary to what is 
usually believed, there is no ambiguity in fact. Thus, any Hermitian construction 
of the quantum Hamiltonian will be necessarily equivalent to that due to any 
other, and consequently nonambiguous. This observation supports the recent 



work of Dutra and Almeida 22 



On the basis of the supersymmetric ideas we have generalized the concept 
of shape invariance to the nonconstant scenario and shown that an appropriate 
choice of the potential and mass variation with the position makes clear the 
link between the effective potentials and the ordering ambiguity. The bound 
state spectra of systems with effective mass are relavant in many areas such as 
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the study of nuclei [28] and metal clusters [ 29[ . However, in other fields such 
as for instance electrònic properties of semiconductors, one is interested in the 
properties of quantum systems with nonconstant mass endowed with continuum 
spectra [1,10,11]. The supersymmetric formalism developed here can also be 
useful tool for the treatment of such kinds of problems and it may be possible to 
generalize our results to such problems. Such an investigation will be deferred to 
a later publication 

As a final remark, an exponentially changing mass and potential have been 
considered in this paper, and it has been shown that the energy levels could 
be redefined in such a way that the ordering ambiguity disappears. However, 
we should note that the exact solvability does depend upon both the form of 
the potential and the change in mass with the position. For instance, as was 
discussed in 0], this cannot be done for the consideration of a quadratically 
growing mass in a singular potential field. However, we believe that in addition 
to its practical applications, the study of quantum mechanical systems with a 
position-dependent mass within the frame of the present technique will raise 
many interesting conceptual problems of a fundemantal nature. In particular the 



method should find wide applications in the study of quasi-exactly solvable |30 



and conditionally-exactly solvable ]3Ï|] systems with nonconstant masses. Along 
this line the works are in progress. 
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Table 1: The special cases of the effective Hamiltonian in Eq. 2. 



Hamiltonian 


a 


a 


P 


7 


Ref. [12] 








-1 





Ref. [13,14] 





-1 








Ref. [15] 





-1/2 





-1/2 


Ref. [16] 








-1/2 


-1/2 


Ref. [21] 
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